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We report an imbibition experiment in 2D random porous 
media in which height - height correlation function grows with 
a nonuniversal exponent. We find the exponent to depend 
on evaporation. A cellular automaton model for imbibition 
based on first principles is presented. A numerical study of 
the model gives results which are consistent with the exper- 
iment. The interface is shown to exhibit self-organised criti- 
cality (SOC). 
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I. INTRODUCTION 

Interface growth phenomena has been studied in a va- 
riety of systems [1-3]. Recently viscous imbibition ex- 
periments with water displacing air in a Hele-Shaw cell 
filled with glass beads opened up a new area to study 
interface growth [4]. Of late, another kind of imbibition 
experiments was studied to understand interface growth 
[5-7]. In these experiments a paper is fixed with the bot- 
tom end dipped into a suspension. The suspension is 
imbibed into the pores of the paper by capillary action. 
The fluid rises through the pores carrying the suspended 
particles with it. Also, as the fluid rises, it continuously 
gets evaporated from the surface of the paper. An inter- 
face is formed by the wet front which rises steadily. There 
is some amount of randomness present in the paper due 
to the blocked pores. This randomness and evaporation 
of the fluid try to pin this wetting front. The motion 
of the wet front is impeded not only by the evaporation 
rate but also by the concentration of the suspension, the 
ratio of the size of the suspended particles to the pore 
size and viscosity of the fluid. The front stops moving 
when the fluid at the boundary has completely evapo- 
rated. These experiments provide a simple laboratory 
system which mimics diffusion in random porous media. 
These are of importance in chromatography. They could 
also provide good "table top " systems to study pattern 
formation [8,9]. 

There are some qualitative differences between this in- 
terface formation in imbibition and other growth models. 
The roughness of the interface in imbibition is solely due 
to the disorder in the paper. If there is no blocking of the 
pores the interface will be smooth at any rate of evapo- 
ration. However, the disorder experienced by the fluid is 
indirectly affected by the evaporation. Remember there 
could be many paths to reach a particular cell. When 
the evaporation is small the fluid can take even longer 



paths to reach a cell. This limits the effectiveness of the 
randomness of the medium in pinning the interface. For 
the same amount of randomness the number of paths will 
be lesser when the evaporation is higher, because higher 
evaporation eliminates longer paths. Thus randomness 
asserts more when evaporation is higher. It should be 
mentioned that this imbibition case differs from the de- 
position models in three ways (i) there is a time scale 
set by the evaporation (ii) this evaporation is from the 
bulk not from the interface alone. Thus evaporation cou- 
ples the interface to the bulk and (iii) the randomness in 
imbibition is quenched and not fluctuating in time. 

It is known that the wet front describes a self-afiine 
fractal roughness [1]. We fix the geometry of the system 
under consideration as follows. Assume that the inter- 
face to be oriented in the x direction. h{x) is the height 
variable along x. The roughness of the interface can then 
be characterised by the exponent a defined as 

W{1) ~ r (1) 

where, W{1) is the correlation of the fluctuations in the 
height at two points x and x + I after the front stops 
moving and is given by 



W{1) = {{h{x)-h{x + l)f) 



1/2 



Till now most of the understanding of the interface 
growth phenomeona in imbibition comes through com- 
puter simulations. So far no simulations in imbibition 
have been done taking care of evaporation, size of the 
particles and concentration of the suspension explicitly. 
In the treatment by Amaral et al [6] evaporation was in- 
corporated phenomenologicaly by a steady increase (Ap) 
at each time step in the probability p of blocking the 
pores. This increase drives the system into a percolation 
threshold pc wherein they obtained a connected cluster of 
blocked pores. This treatment assumes a linear increase 
of p with time. In an alternative description one can as- 
sociate a dimensionless driving force f = — Fc)/Fc 
with the interface, with F decreasing to Fc being equiv- 
alent to p — > pc. That is the change of p in this model 
is in effect a linear decrease of this driving force / with 
time, driving it to Fc- By incorporating evaporation ex- 
plicitly we show that a nontrivial nonlinear dependence 
of driving force will be more appropriate. 

It is believed that in general interface growth can be 
described by a nonlinear growth equation [10,11] 

dth{x) = F + iyV^h{x) + X{Vh{x)Y + ({x, h) (3) 

where, C,{x, h) is quenched noise and F is a driving force. 



1 



The above equation describes two universality classes 
depending on the value of A as F — > F^. In the case of 
Directed Percolation Depinning (DPD) A — > co while in 
quenched Edward - Wilkinson (QEW) universality class 
A — > as F — > Fc. The model treated by Amaral et. 
al. [6] falls into the DPD universality class. However, 
we find imbibition to be different from these two classes 
by the presence of a driving force which tunes itself close 
to Fc leading to self-organised crticality (SOC). 

In general SOC is characterised by the emergence of 
scale invariance in an extended nonequilibrium system 
[12]. Even though there is a substantial amount of work 
in this field till now there is no general aggrement on the 
origin and characterisation of SOC in different systems. 
In fact conservation laws were believed to be necessary 
for the appearance of SOC [13]. However there are also 
models which have nonconservation but still shows SOC 
[14,15]. The model presented here belongs to this class 
of nonconserved systems exhibiting SOC. 

In this paper we present a detailed analysis of the re- 
sults communicated through a short paper [7]. The pa- 
per is organised as follows. In section 2 we describe an 
imbibition experiment at two different experimental con- 
ditions. In section 3 we present two microscopic models 
to study imbibition. A numerical study of the models 
and comparison of the results with that of experiments 
is given in section 4. We point out an interesting connec- 
tion beteween imbibition and other nonconserved models 
and finally in section 5 we give our conclusions. 



II. EXPERIMENT 



The experiments are carried out using Whatman No:l 
filter paper as the porous medium and ink as the suspen- 
sion. The evaporation rate was varied by changing the 
room humidity and temperature. The ink rises through 
the paper and stops at a particular height. The darkening 
of the interface is indicative of this stoppage of growth. 
This darkening is due to the following reason. As the ink 
rises the fluid evaporates. As a result the suspended par- 
ticles carried by the fluid get deposited at the interface 
continuously causing a darkening of the boundary. The 
paper was taken out from the suspension and dried when 
this darkening occurs. The interface was then digitized 
using a CCD camera and a frame grabber with a resolu- 
tion of 260 pixels per inch (see figure 1). Since changing 
the room condition is found to alter the roughness of the 
interface even after it is dried, the humidity and tem- 
perature of the room was kept constant till the image is 
digitized. Figure 2 shows the behavoir of W{1) against 
I for two different values of evaporation rate. The data 
are averaged over 10 experiments. 
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FIG. 1. (a) A digitized image of the interface from exper- 
iment, (b) A typical result of the simulation obtained using 
model I 



FIG. 2. The experimental values of height - height correla- 
tion function W(l) plotted against the distance of separation 1 
after the interface stopped growing. Points marked with * fall 
onto a curve with exponent .45 ± .004. For lower evaporation 
the points (marked o ) correspond to an exponent .67 ± .004 

We find the value of a to be .45 ± .004 and .67 ± .004 
for higher (humidity 50.0% and temperature 21.0°C)and 
lower evaporation (humidity 60.0% and temperature 
23.0°C)respectively. However, we do not have a quan- 
titative measure of the evaporation. Thus a is very dif- 
ferent for the two cases. A dependence of a on the ex- 
ternal parameter is also exhibited by DPD models when 
the growth is stopped before it reaches p^. There the 
exponent crosses over to the KPZ value. The change in 
a shown by the experiments above is different from this 
cross over since we see values of a below KPZ [a = .5) 
and above DPD (a = .63). This shows that unlike the 
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results of reference [6] a is not universal and thus cannot 
be described by DPD model. 

III. THE MODELS 

To understand the dependence of the various parame- 
ters affecting the growth of interface we present two mod- 
els for studying the imbibition of a suspension into 2D 
random porous media. Before discussing the models used 
in the simulation let us try to understand the problem 
from a microscopic point of view. In the experiments 
described above paper was used as the random medium. 
At a microscopic level one can regard the paper as a ran- 
domly disordered medium [6] with a quenched random- 
ness given by a fixed probability p for the pores to be 
blocked. So, the interface growth phenomena is nothing 
but the propagation of fluid particles through this disor- 
dered medium. The wetting front of the fluid particles 
propagates due to the capillary forces. The disorder in 
the medium and the evaporation tries to pin this growth. 
It should be noted that for p < Pc this pinning is pos- 
sible only when the evaporation is present. Evaporation 
constantly decreases the number of fluid particles in the 
wetting front. This makes it more difficult for the fluid 
to overcome the obstacles. The front stops when the 
number of fluid particles there goes to zero. The average 
number of fluid particles in the front can be taken as a 
reduced driving force f = — Fc)/Fc. The front stops 
moving when / — > at a critical height. It is true that 
smaller the evaporation rate, larger will be this critical 
height. 

A. Model I : Cellular Automaton 

In this model the porous medium is considered as a 
square lattice with disorder being incorporated by block- 
ing some cells randomly with a constant probability 
p < Pc (see flgure 3). The maximum capacity of each 
cell is flxed to Nq number of fluid particles i.e. the cells 
are of flxed volume and the fluid is incompressible. At 
every time step, evaporation was explicitly modelled by 
the loss of certain number of particles n in the trans- 
fer. At time t=0, at the bottom edge of the lattice a 
horizontal line of wet cells with Nq particles is created. 
At t+1 the particles are imbibed into all unblocked cells 
which are nearest neighbors to the wet region. If N is 
the number of particles in a wet cell then it transfers 
N — n particles to all its unblocked nearest neighbours 
I.e. each of the neighbours get N — n particles, n being 
the loss due to evaporation. This updating is done par- 
allely i.e. all cells which are nearest neighbour to the 
front are updated simultaneously. When a cell transfers 
to its nearest neighbours the number of particles it con- 
tains remains the same due to the source below. If a cell 



has more than one wet nearest neighbour it gets particles 
from all of them subject to a maximum number Nq. We 
also apply the rule that every cell blocked or unblocked 
below a new wet cell become wet as well [6]. This is to 
avoid the presence of overhangs and islands. We use pe- 
riodic boundary condition in x by identifying the cells at 
the edges of the lattice. 
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FIG. 3. Example of the multiple connectivity of the model 
for a 6X6 lattice. The blocked cells are shown black. At time 
all the cells corresponding to i=l and j = l,6 are filled with 
No particles. At each time step the particles are transfered 
to nearest neighbours. Note that the cell (4,4) gets particles 
from both (4,3) and (4,5). Also when a cell (is,j) is wet all 
the other cells (i,j) with i < is are wet as well. 

The concentration of the suspension and the pore size 
are incorporated in the model in the following way. For 
a given concentration the ratio of the number of fluid 
particles to the number of suspended particles in a cell is 
flxed. Since Nq is the maximum allowed number of fluid 
particles in a cell and a cell is fllled with both the fluid 
and suspended particles a change in concentration of the 
suspension effectively changes the number Nq. Thus a 
change in concentration of the suspension can be mod- 
elled by changing the maximum allowed number of fluid 
particles in a cell Nq. 

In this model we can apply an external bias which 
could be present due to gravity or anisotropy in the 
medium. This bias was incorporated in the model by 
introducing a difference in the number of particles trans- 
fered to the vertical and horizontal neighbours of a given 
cell. 

There are two independent parameters in the model. 
They are (i) the probability p and (ii) the ratio 77 = 
n/No. 

The height-height correlation function W[l) in this 
model exhibits a power law with a nonuniversal exponent 
a which changes continuously with 77 [7]. The model is 
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different from that of reference [6], in the sense that here 
the effect of evaporation is incorporated in an explicit 
way. Unlike the directed percolation models where the 
interface is pinned by the connected cluster of blocks here 
it stops only when the wetting front runs out of fluid due 
to evaporation. When this happens a connected cluster 
of dry cells (which could be blocked or unblocked) forms 
at the boundary. It should be noted that the running 
out of fluid is qualitatively different from stopping the 
DPD growth before it reaches p^. This is clear from the 
fact that for all values of evaporation the interface stops 
forming a connected cluster of dry cells. The average 
number of particles in the wetting front can be consid- 
ered as the reduced driving force /. Thus evaporation 
continuously changes this driving force. Depending on 
the nature of the change in driving force as a function of 
time, the growth behaviour of the interface also changes 
qualitatively. The model is also different from the Eden 
growth [1] in the sense that all the sites in the boundary 
moves at the same time i.e. here the growth process by 
itself does not cause interface roughness. 

The model can be used to study, for example, thin layer 
chromatography wherein a mixture of different chemical 
compounds are made to diffuse through a porous medium 
resulting in their separation. 

B. Model II : Random Walk 



The function < PiPi+i >j varies as P. The exponent 
a in this model is a continuous function of the death rate 
d. We find that there is a mapping between model-I and 
model-II with d O p and 1 — 6 O The equivalent of 
reduced driving force / here is the number of live cells at 
the boundary. All the qualitative features of model-I can 
also be obtained from model-II. In the discussion below 
we consider in detail model-I only. Another interesting 
feature of imbibition brought out by this model is it's 
close similarity with reaction diffusion systems which has 
been studied in great detail by many authors [16]. 

IV. RESULTS AND DISCUSSION 

In figure 4 we show the behaviour of W{1) as a func- 
tion of I for different values of evaporation for a fixed 
value of p and lattice size L. The simulations show the 
existence of a crossover length such that the height - 
height correlation function VF ~ P for Z <C Ix- Whereas 
for I ^ lx,W saturates to a constant value Wsat- The 
exponent a is a function of evaporation 77 while the value 
of Wsat and depend on both the system size L and 
evaporation 77. The dependence of a and Wsat on evapo- 
ration is consistent with that found in the experiment as 
can be seen by comparing figures 2 and 4. 



The qualitative behaviour of the interface is not very 
sensitive to the microscopic details of model I. This can 
be demonstrated by the fact that a random walk model 
with a finite birth - death probability captures all the 
features of imbibition qualitatively. The model is as fol- 
lows. Consider a variable pi. This can be thought of as 
the occupation number or the height of each column i in 
a square lattice. To start with we set all pi = 0. The cells 
at the tip of the columns (boundary cells) are considered 
as random walkers. To make exact correspondance with 
model-I we update all the columns simultaneously. The 
boundary cells either die or diffuse with fixed rates d and 
1 — d respectively. 

If diffusion is chosen, then one of it's nearest neighbour 
is picked as the target with equal probability. If a cell up 
or down the boundary cell of column i is chosen, then pi 
is changed by ±1. However, if the picked cell is one at left 
or right then two events are possible. If the neighbour 
is a dead one then it is regenerated with a probability b. 
In the case of alive neighbour the value of p for both the 
columns are set equal to the larger one. 

If death was chosen instead of diffusion then the value 
of Pi corresponding to that column is kept the same. This 
column is then not considered for updating untill regen- 
erated by one of its live neighbours. The growth process 
stops when all the boundary cells are dead. 
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FIG. 4. The simulated height-height correlation function 
W^l) plotted against the distance of separation I after the 
columns stopped growing. The length is measured in units 
of lattice parameter. The parameters are p = .45, 7; = .12 
to .148 in equal intervals of .004. The saturation value of 
W^l) decreases monotonically with 7;. The simulations are 
done with a lattice of length 6000 applying periodic boundary 
condition. 

The exponent /3 defined through the time correlation 
function C(t) as 
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C{r) = {{h{t + r)-h{t)Y) 



2\l/2 



(4) 



is shown in figure 5 for various values of 77. Note that like 
a, (3 also exhibits a dependence on evaporation. 




FIG. 5. The simulated time correlation function of height 
for values of 7; = .120, .133, .140, .148 and other parameters 
same as figure(4) 

As mentioned earlier a crossover of a from .63 to 0.5 
is seen in DPD model [5]. However, the dependence of 
the exponents on evaporation seen in the present model 
is different from this crossover. This is supported by the 
following arguments, (i) In both simulation and experi- 
ment we find the exponent a to vary from a value below 
0.5 (KPZ) to one above 0.63 (DP), (ii) The dynamic ex- 
ponent (3 obtained from the simulation is not the same 
as the KPZ value in the region where a = 0.5 (iii) Un- 
like the crossover region observed in DPD models here 
the growth is stopped only when F = Fc (iv) We find 
a to vary continuously with 77, i.e. we do not find any 
wide enough region of 77 where scaling with a single value 
a can be obtained and finally (v) KPZ is expected when 
the interface is moving with a constant velocity, like driv- 
ing DP model with p < pc. But in the present model even 
though p < Pc evaporation continuously changes the ve- 
locity. These points show that the underlying physics 
here is different from that of the DPD and KPZ models. 

To get an insight into the dependence of a and Wsat on 
the parameters mentioned above we plot in figure 6 the 
change in reduced driving force / as a function of time. 
It is shown that for the range of 77 considered / behaves 
in the following form 

/(i) ~exp[^/(i + i,)]. (5) 

Here A = A[r]), is a function of evaporation which sets 
the initial time scale. This does not affect the long time 
behaviour, tc is a constant for a given value of p. Chang- 
ing p does not change the behaviour of / qualitatively but 



will change the values of ^ and tc. This nonlinear depen- 
dence of / on time explains why the interface behaviour 
is different from that of DPD model in ref [6]. 
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FIG. 6. The nature of the reduced driving force / for 
71 = .136 and p = .45. The continuous line is the fit to equa- 
tion (5) with tc = 10.0 and A = 29.0 

The driving force / saturates to a value /j at i ~ 2tc. 
The value of /j as a function of 77 is given in figure 7 
at t = 2tc. We would like to emphasize the point that 
changing 77 do not change tc but merely changes the value 
of fc at which the system is driven and A which sets the 
transient time scale. We find that there are two critical 
values of 77, 771 and 772 at which the slope of the /s(77) curve 
changes drastically. The interface behaves differently in 
this three regions. 
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FIG. 7. The reduced driving force fs = /(2ic) showing two 
critical values of 7; for p = .45 

It is easy to understand the reason for the existence of 
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two critical values of 77. In the model the particles have 
more than one path to reach a particular cell. The effec- 
tive number of paths available to reach a cell decreases 
with evaporation since increased evaporation suppresses 
the longer paths. For ?7 < ?7i the number of fluid par- 
ticles that a cell loses through evaporation is more than 
compensated by the inflow because of the many paths 
available. Hence this regime becomes super diffusive. We 
see from figure 8 that the value of a in this region to be 
greater than .5. As the number of paths become less 
the particles get stuck at obstacles for a longer time in- 
ducing a transition into a normal diffusive regime with 
exponents close to .5. In this region the reduced driving 
force / saturates to a non-zero value (/j ~ 0+) at i ~ 2tc. 
This driving of the interface close to / ~ 0+ gives rise to 
some interesting behaviour. We will discuss this region 
in detail separately. On further increase of evaporation 
the driving force / goes to zero at i = 2tc. The system in 
this region becomes subdiffusive with exponents falling 
much below .5. 
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FIG. 8. The dependence of the simulated exponent a on 7; 
at i = 4ic for p = .45 



A. Scale invariant region (t/i < t; < 7/2) 

This is the most interesting region. Here W{1) shows 
a dependence on the system size L. This is depicted in 
figure 9. We find that the cutoff width Wsat depends on 
the system size as 

W,,t ~ L\ (6) 

with the width function W{1, L) satisfying a finite size 
scaling form 



W{l,L) = rg{l/L''). 



(7) 



Where 7 = .34 and v = .68 for p = .2 and 77 = .293 
The data collapse obtained with this scaling is shown in 



figure 10. As implied by equations (1) and (7) we find the 
exponents 7, ly and a satisfies the scaling relation [17] 
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FIG. 9. The simulated height-height correlation function 
W^l) plotted against the distance of separation I after the 
columns stopped growing for system sizes L = 125, 250, 500, 
750, 1000, 1250, 1500, 1750,2000, p = .2, 7; = .293 
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FIG. 10. The data collapse with the finite size scaling (see 
equation (7)) for 7 = .34 and u = .68. The other parameters 
are same as in figure 9 
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since for the given value of p and r] a = .52. This system 
size scaling establishes that there is no single character- 
istic length scale in the problem. From our simulations 
for the same value of p we find another relation 



1=1 
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However the scaling relation 7 + = /x obtained from 
equations (7) and (9) is not obeyed here pointing to mul- 
tifractal behaviour. This implies that the dominant na- 
ture of the function W[l, L) is the power law relation 
given by equation (1) [17]. 

As mentioned earlier the interface stops moving when 
all the boundary cells dry up. Then a connected cluster 
of dry cells is formed at the boundary. We show in figure 
11 that the length distribution of this connected clusters 
approches a power law distribution, 

p(0 ~ (10) 

as i — > CO. This power law dependence is typical of 
scale invariant systems. 
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log(e) 

FIG. 11. Distribution of the length of connected cluster 
of dry cells calculated at i = ic to i = Sic- Note that the 
distribution approches to a power law given by P(^) ~ ■ 

B. self-organised criticality (SOC) 

To establish the SOC nature of the interface in the 
above models it is necessary to obtain the length and 
time distribution of avalanches [18]. We get this by the 
following method. When the interface stops moving a 
connected cluster of dry cells is formed at the boundary. 
We then fill one of these dry cells with Nq, number of par- 
ticles. This results in further movement of the interface. 
The spatial and temporal distribution of these avalanches 
is shown in figure 12. These distributions are found to 
obey the power laws 

P{1) ~ l-"- (11) 
P{t) ~ i"*. (12) 
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FIG. 12. Spatial and temporal distribution functions of the 
avalanches shown for the parameters same as in figure 6 

The exponent 6 has a universal value equal to 1.0 for 
rii < T] < ri2. It has been argued recently by Boer et 
al. that the mere existence of power law correlations is 
not sufficient to establish SOC [18]. They take the non- 
normalisability of the temporal distribution function to 
be the criterion of SOC. The corresponding argument 
here is the following. For an infinite lattice the finiteness 
of the integral P(i)di means that points infinitely far 
from the origin of avalanche do not participate in the 
avalanche process. If the interface is critical it is robust 
against any disturbance, i.e. an avalanche should leave 
the interface unchanged. This is possible only if all the 
points on the interface has a nonzero probability of par- 
ticipating in the avalanche indicating P{t) as i — > co to 
be significant. This results in the divergence of the in- 
egral. The value of 5 = 1.0 in our simulation satisfies this 
criterion establishing the SOC of the interface. 

There exists a simple relation between the exponents 
K and a given by k = 1.0 + a . This can be seen from 
the following. The total area of an avalanche of size I 
is J^l°'dl ~ 1^'^°'. Since there is a uniform probability 
of blocking the probability of the avalanche of size I is 
;-(i+a)_ The same argument applies for the distribution 
P{C) shown in figure (11). 

To understand this critical behaviour of the interface 
we look at the nature of the driving force. Remember, 
(see figure 6) the time dependence of the driving force is 
of the form given in equation(5). This means that the 
interface is driven by a force F ~ Fc for a long time. It 
should be noted that the tuning of the driving force to Fc 
is done by the system internally. This is a characterestic 
feature of self organising systems, wherein the system is 
driven in a region close to criticality by itself [19]. This 
shows that in this region of 77 the interface exhibits self 
organisation. 

The behaviour of correlation function W[l) on evapo- 
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ration is similar to the dependence of distribution func- 
tions on nonconservation in cellular automaton in earth- 
quake models [14]. In fact the dependence of the expo- 
nent a on evaporation is strikingly similar to that of the 
exponent B on nonconservation parameter in reference 
[14]. This indicates that imbibition is better described 
by cellular automaton models than directed percolation. 
The main points which distinguishes the present model 
from that of DP are the following, (i) Unlike the present 
model DP model has universal exponents, (ii) the DP 
model is driven to pc in a linear fashion whereas the 
present model has a nonlinear driving force leading to 
a terminal driving of the interface and most importantly 
(ii) we find the avalanches to show power law distribution 
while the avalanches in the DP model shows no interest- 
ing scaling behaviour and has an exponential distribution 
[20]. This power law distribution shows that evaporation 
does not set any length scale contrarary to the DP model 
wherein there is a length scale set by the value of Ap. 



C. Multifractality 



of height of the columns from the mean h. We define a 
normalised probability distribution 



ET=-M m 

Then the singularity spectrum f[6) is defined as 



(14) 



m 



ln(2M 



1 ^ 

^ E mi)invr(i) (15) 



j = -M 



where, 



and 



-1 ^ 

e = — V invr(i) 

ln(2M + 1) 



j=-M 



To check for multiscaling of the interface we calculate 
from our simulation the mth order height-height correla- 
tion function defined by [21] 



Wrn{l) = {\h{x + I) - h{x)r) ^ F 



fl3) 



The value of the exponent for various values of m 
is shown in figure 13. As is well known for multifractal 
surfaces we find to change continuously with m. 
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FIG. 14. /(5) versus 6 spectrum for the same parameters 
as figure 13 

In figure(14) we depict the smooth dependence of f[6) 
on 6. Remember for ordinary fractals this curve degen- 
erates into a single point. 

It is to be noted that the results are presented for a 
particular choice of parameters in the model. A change 
in the parameters does not alter the gross features. 



m 

FIG. 13. The exponent a-m plotted against m. Parameters 
are p = .45 and ri = .136 

We also carry out a multifractal analysis using the 
height distribution N[j) [22,23]. Here j is the deviation 



V. CONCLUSION 

The main points of this paper are the following. 

(i) We have shown from experiments conducted at two 
different conditions that the exponent a in imbibition is 
not universal. 
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(ii) A model for imbibition (model I) which includes the 
effect of evaporation through loss of particles on transfer 
is presented. 

(iii) The model shows that the static height - height 
correlation exponent a depends on the evaporation which 
is consistent with our experiments. 

(iv) A random walk model (model II) which captures 
all the qualitative features of imbibition is presented to 
show the independence of the general results on micro- 
scopic details. 

(v) The interface is shown to exhibit self-organised crit- 
icality for a region of evaporation. The dependence of the 
exponents on evaporation is simiar to that of exponents 
in nonconserved cellular automaton indicating that im- 
bibition falls into this class. 

(vi) It is shown that the driving force tunes itself close 
to the critical value for all values of evaporation 771 < 
?7 < ?72. This slow driving of the system close to pinning 
can be understood as the reason for the existence of self- 
organisation [19]. 

(vii) The interface also shows multifractality. 

To conclude, imbibition experiments show that the 
roughness exponents depend on evaporation. A first prin- 
ciples cellular automaton model for imbibition gives re- 
sults in agreement with experiments. Also the model 
exhibits properties much richer than the already known 
models of imbibition. 
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